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Abstract 

Six-dimensional Einstein-Gauss-Bonnet gravity (with a linear Gauss-Bonnet term) is inves- 
tigated. This theory is inspired by basic features of results coming from string and M-theory. 
Dynamical compactification is carried out and it is seen that a four- dimensional accelerating 
FRW universe is recovered, when the two-dimensional internal space radius shrinks. A non- 
perturbative structure of the corresponding theory is identified which has either three or one 
stable fixed points, depending on the Gauss-Bonnet coupling being positive or negative. A much 
richer structure than in the case of the perturbative regime of the dynamical compactification 
recently studied by Andrew, Bolen, and Middleton is exhibited. 
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1. Introduction. To realize that the expansion of the universe is accelerating was one of the 
most important scientific discoveries of the last century. There are several alternative explana- 
tions of this remarkable fact (what already means that it is not well understood yet). A quite 
appealing possibility for the gravitational origin of the dark energy responsible for this acceler- 
ated expansion is the modification of General Relativity (GR) or the corresponding Einsteinian 
gravity (see [IJ for a review of these approaches). Such a well-established and successful theory 
cannot be modified without a very good reason, and much less in an arbitrary way. But observe 
that there is no compelhng reason why standard GR should be trusted at cosmological scales. 
For a rather minimal modification of the same, one assumes that the gravitational action might 
contain some additional terms which would start to grow slowly with decreasing curvature. A 
particularly interesting formulation is obtained by using well-grounded geometrical arguments, 
specifically when the modified gravity action is endowed with a function of the Gauss-Bonnet 
(GB) topological invariant, G, as it was suggested in [2]. 

It must be noted that different types of dark energy may actually show up in different ways, 
at large distances. Cold dark matter is known to be localized near galaxy clusters but, quite on 
the contrary, dark energy distributes uniformly in the universe. The reason for that behavior 
could be explained by a difference in the equation of state parameter w = p/ p. Moreover, the 
effect of gravity on the cosmological fluid turns out to depend on w and it so happens that, even 
when — 1 < w < 0, gravity can act sometimes as a repulsive force. The effect of gravity on matter 
with — 1 < w < can be shown to be opposite to that on usual matter, which becomes dense 
near a star, while matter with — 1 < w < becomes less dense when approaching a star [2]. 
Dark energy contributes uniformly throughout the universe, which would be indeed consistent, 
since the equation of state parameter of dark energy is almost —1. If dark energy is of phantom 
nature {w < —1), its density becomes large near the cluster but if dark energy is of quintessence 
type (—1 < w < —1/3), then its density becomes smaller. 

Another very important argument in favor of the GB modified theory is the fact that it can be 
seen as being inspired by string and/or M-theory [3]. In fact, specific models for string-inspired 
scalar GB gravity, considered as possible forms of dark energy, have been discussed in [1] and 
[5]. It was subsequently shown in [6], that scalar GB gravity can actually be represented as a 
modified GB gravity without scalars, and specifically that it can be equivalent to an ideal fluid 
with an homogeneous equation of state [7]. 

In the present paper we will study the case of 6-dimensional Einstein GB gravity, with a 
linear GB term. As already advanced, this situation is quite interesting, since the theory that we 
will thereby obtain can be shown to be inspired by what is derived from string/M-theory, after 
some specific compactification is carried out, when the scale factor of the 2-dimensional internal 
space goes to zero. In the course of our study we will discover a non-perturbative structure of 
the corresponding theory with either three (for positive GB coupling, e) or one (for negative GB 
coupling) stable fixed points. This will exhibit a much richer structure than the one that follows 
from the case of the perturbative regime of the dynamical compactification recently studied by 
Andrew et al. [8] 

2. Dynamical compactification of 6-dimensional Einstein-Gauss-Bonnet gravity. We 

shall start from the following, string-inspired action in six dimensions 



where e is a constant, and the metric is the product of the usual metric corresponding to the 
4-dimensional FRW universe and a 2-dimensional surface, namely 





ds^ = -de + a\t)[{dx^f + {dx'f + [dx^'f] + h''{t)[{dx^f + {dx^'fl 



(2) 
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the scalar curvature is 

^ 6(2^ 12ab 26^ 6a 46 
ab b'^ a 

while the four-dimensional and topologically invariant Gauss-Bonnet Lagrangian, Lgb-, has the 
form 

^ _ ASa% ^ 72a%'^ ^ 24a^'d ^ 96adb ^ 2A'db'^ _^ A9,d% ^ 48abb 
a% a?b'^ a?b ah^ a?b ah^ ' 

or, equivalently, 

24 . . . . ..... 

Lgb = -^{2(x%b + 3aa^6^ + d^afe^ + Aaixdbb + a^dfe^ + 2a(x%b + 20^ abb). (5) 

Note that, unlike the six-dimensional compatification case of Elizalde et al. [9J, the theory under 
consideration is not multiplicatively renormalizable. The corresponding equations of motion are 
obtained by variation of the action with respect to a and b, what yields 

d^6^ + Aadbb + a%'^ + 2ad6^ + 2o?bb + I2ea%'^ + IQeadbb + 8ead6^ 

+Sta?bb + IQtaabb = 0, 
3ad^6 + 3a^d6 + 3a^d6 + a% + 12ed^6 + Uea^db + 2Aeaddb + 12ead^6 = 0. (6) 

These equations can be easily rewritten in terms of the Hubble rates H = a' /a and h = b'/b, 
namely 

3/i^ + AhH + SH"^ + 2h + 2H+ IGeh^H + 28eh^H'^ + IGehH^ 

+16ehhH + 8ehH^ + 8eh^H + 16ehHH = 0, 
+ 3hH + 6H^ + h + 3H + Ueh'^H'^ + 36ehH^ + 12eH^ 

+12ehH'^ + 2AehHH + 12eH^H = 0. (7) 

In addition, variation over the metric in the above expressions gives the constraint equation 

h'^ + QhH + 3H^ + 3Qeh'^H^ + 2AehH^ = (8) 

This equation helps to exclude h and h' from Eqs. ([7]). As a result, one gets an equation for H 
only: 

H' = 3H'^ X (9) 
VQ + AG + e{-22V6 + Q4G)H'^-24e\9VQ - 52G)H%9Qe^{17VQ + 12G)H^-8064VQe^H^ 
VQ - 12e{VQ + IQG)H^ + 72e2(3V6 - 32G)m - 288Qy/Qe^H^ + 3llQAy/Qe^H^ ' 

where 



G = Vl - 6eH^ + 2Ae'^m. (10) 
One can check that this last equation obeys the fundamental relation (for e > 0): 

^ - (iJ2_^2)(^2_,2)/(^)' (11) 
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where g and r are constants, and the function f{H) < 0. We start now with its numerical 
analysis. 

3. First case: e > 0. It leads to the following values of the constants, corresponding to constant 
curvatures in the 4-dimensional and 2-dimensional spaces: 

po ~ 0, 7501/v^, go ~ 0,4842/v^, ro ^ 0, 1023/v^. (12) 

Note that po is a stable fixed point while go and vq are singular points. It turns out then, 
that the initial values of H{t) can be classified as belonging to four different regions, which are 
delimited by these values of po, go and tq. The behavior of H(t) in each of these regions can differ 
considerably, from one to the other. We will now consider the different possibilities in detail. 
The four different cases corresponding to e > are as follows. 



1. Case < H{0) < tq, then H' > 0, when H — * go, and in the limit it turns out that H' +oo. 
Thus, at if = go a singularity develops. Figs, la, lb (for a classification of the different types of 
future singularities, see e.g. [TU]): 
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Fig. Ic Fig. Id 

We see that in this case the resulting FRW universe expands with acceleration while the com- 
pactification radius of the extra-dimensional space decreases with time. This indicates that the 
higher-dimensional GB term can indeed play the role of dark energy in this universe. As is the 
cases in the majority of the dark energy models currently available, a future singularity occurs. 

2. Case tq < H{0) < go, then H' < 0, when H — »• ro, and in the limit one has that H' — > —oo. 
Again, H = tq corresponds to a singularity, and the curves are, respectively (Figs. 2a, 2b): 



4 



. 4 
K 0.3 




0.5 1 1.5 2 0.1 0.2 0.3 0.4 

Fig. 2a ^ Fig. 2b " 

As in the case before, using Eq. ([8]) the behavior of h{t) and h'{h) can be obtained (Figs. 2c, 2d) 
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Fig. 2c Fig. 2d 

We see that in this case the size of the FRW universe decreases, while the 2-dimensional inter- 
nal space scale factor may increase. This case cannot thus naturally describe the dark energy 
universe. 

3. Case go < H{0) < po, then H' > 0, when H Pq, and then H tends to a constant, after an 
oscillatory regime (Figs. 3a, 3b). 




Fig. 3a Fig. 3b 

Furthermore, the scale factor of the internal space can actually decrease, in fact (Figs. 3c, 3d) 
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Fig. 3c Fig. 3d 

Hence, the exact solution obtained at H = po turns out to be a stationary stable point. Precise 
analysis corresponding to specific values of the parameters can be further carried out in a rather 
simple way. 

4. Case < H{0), then if' < as if — > po, and the solution in this case is oscillatory (Figs. 4a, 
4b). 




Fig. 4c Fig. 4d 

In this situation the FRW scale factor tends to a stationary point with decreasing H. The dark 
energy universe can correspond, in this case, to one of the branches of the oscillatory universe: 
expansion is turned into contraction, and vice- versa, with the repeated oscillations. Such universe 
would explicitly correspond to the general case described in [TT]. The analytic solution for this 
case can be obtained. 

Thus, it turns that in the FRW universe the scale factor behavior is much less complicated 



6 



for these models than in the case of the absence of the GB term. This is certainly rewarding. In 
particular, for instance, in order to keep close to the Einstein regime, the condition of proportion- 
ality (perturbative regime) of a{t) and b{t) was used for the study of dynamical compactification 
in the Einstein-GB theory in Ref. |8J. For comparison, one can easily check that the situation 
there gets quite complicated and that different regimes for the Hubble rates appear, everything 
being much more simple and natural in our model above. 

4. Second case: e < 0. In this case there is one stationary stable point H = h = ±1/ \/6|e|. 
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Fig. 5b 
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Fig. 5c 
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Notice here that solutions with close initial data tend to the stationary point in the oscillatory 
regime. However, both scale factors have the same sign and no dynamical compactification 
occurs! The universe, as a whole 6-dimensional object, either exponentially expands or shrinks. 
This indicates already that problems should be expected from a negative— e model (indeed, in 
the higher-dimensional black hole case this may lead to negative entropy |il4j). 

5. Third case: e = 0. In this case one does recover (as it should be) an explicit solution. In 
fact, the equations for a, b are 

b'^a'^ + Aaba'b' + a%'^ + 2ab'^a" + 2a%b" 



Saba'"^ + Sa^a'b' + 3a^ba" + a^b" 

3h^ + AhH + SH"^ + 2h' + 2H' 
+ 3hH + QH"^ + h' + 3H' 

H' ={3± 2VQ)H^ 



and, from here. 



From where one gets that 



0, 
0, 



0, 
0. 



(13) 

(14) 
(15) 
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and the solution is given by 



H = 



1 



(16) 



being 

a = 3±2V6. (17) 

Moreover, in terms of the scale factors: 

a = C2[±(at + Ci)]-^/", 

b = C3[±(at + Ci)]^/", (18) 

where 

P = 3±V6. (19) 

6. Discussion. Summing up, we have investigated in this paper explicit non-perturbative dy- 
namical compactification to a 4-dimensional FRW universe starting from a model of 6-dimensional 
Einstein-GB gravity. The number of stationary points obtained depends on the sign of the GB pa- 
rameter, e. They correspond in fact to exact solutions where the curvatures of the 2-dimensional 
and 4-dimensional spaces are constant. We have found a regime where the FRW universe does 
expand with acceleration, at the same time that the scale factor of the internal space goes to zero. 
But the most interesting regime discovered here is the oscillatory case, where a unification of an 
inflationary epoch with a late-time acceleration regime is indeed possible. This was suggested 
in [12] and has been explicitly realized here with concrete solutions and numbers, by combining 
numerical with analytical methods. For lack of space, in this letter we have just provided a brief 
presentation of the results. A more detailed study of the different accelerating universes that 
are derived from the model will be given elsewhere, in particular, details on the sequence of the 
matter dominated phase, the transition to acceleration, and the proper accelerating phases, as 
was done for other dark energy models of modified gravity in [13]. It will be also necessary to 
include in our considerations the effects of scalars (dilaton, moduli) in the low-energy sector of 
the string effective action. 
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